THE WIDTH OF RESONANCES FOR SLOWLY VARYING 
PERTURBATIONS OF ONE-DIMENSIONAL PERIODIC 
SCHRODINGER OPERATORS 

FREDERIC KLOPP AND MAGALI MARX 

Abstract. In this talk, we report on results about the width of the resonances for a 
slowly varying perturbation of a periodic operator. The study takes place in dimension 
one. The perturbation is assumed to be analytic and local in the sense that it tends to a 
constant at +00 and at —00; these constants may differ. Modulo an assumption on the 
relative position of the range of the local perturbation with respect to the spectrum of 
the background periodic operator, we show that the width of the resonances is essentially 
given by a tunneling effect in a suitable phase space. 

Resume. Dans cet expose, nous decrirons le calcul de la largeur des resonances de 
perturbations lentcs d'opcratcurs de Schrodinger pcriodiqucs. Cette etude est uni- 
dimcnsionncllc. Lcs perturbations lentes considerees sont analytiques et locales au sens 
ou elles tendent vers une constante en +00 et en — 00 ; ces deux constantcs peuvent toute- 
fois etre differcntes. Sous des hypotheses adequates sur la position relative de l'image de 
la perturbation locale par rapport au spectre de l'operateur de Schrodinger pcriodiquc, 
nous demontrons que la largeur des resonances est donnee par un effet tunnel dans un 
espace de phase adequat. 



0. Introduction 

The present talk is devoted to the analysis of the family of one-dimensional quasi- 
periodic Schrodinger operators acting on L 2 (R) defined by 

(0.1) H c , £ = + V(x) + W(ex + C). 

We assume that 

(HI): V : M. — > M. is a non constant, locally square integrable, 1-periodic function; 
(H2): e is a small positive number; 
(H3): ( is a real parameter; 

(H4): W is a potential that is real analytic in a conic neighboring the real axis that 



admits a limit at +00 and at —00; the precise assumption is stated in section [TT2 . 



As e is small, the operator ( p.l|) is a slow perturbation of the periodic Schrodinger operator 
(0.2) Ho = -^ + V(x) 
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acting on L 2 ( 

When V = 0, the operator H^ E is independent of ( (up to a unitary equivalence, a trans- 
lation actually) and becomes a semi-classical Schrodinger operator (by a simple rescaling 
of the variable x). The resonances for such operators have been the subject of a vast 
literature over the last 20 years; a detailed review and many references can be found in 
the lecture notes and in the papers . 

The problem of computing the resonances for slowly varying perturbations of a periodic 
Schrodinger operator has been much less studied. The papers H | deal with the multi- 
dimensional case; asymptotic formulas are obtained for the real part of the resonances. 
This is done using trace formulas and the imaginary parts of the resonances have not been 
computed. The papers |], |2| are devoted to the study of Stark- Wannier resonances in a 
small electric field. 

It is well known that the eigenvalues and resonances for H^ E near an energy E 
depends very much on the region of energy one is studying. Les us assume W tends to 
at both ends of the real axis. Then, the absolutely continuous spectrum of H^ E is the 
spectrum of Hq. In this case, we compute the resonances of H^ £ near energies inside 
the spectrum of Hq. Actually, for W, we also consider the case when the limits of W 
at +oo and — oo differ. We obtain that the real parts of the resonances are given by a 
quantization condition interpreted naturally in the adiabatic phase space (see section |iT6|) . 
The imaginary parts of the resonances are then given by an exponentially small tunneling 
coefficient. 

Though the results depend crucially on the case dealt with, from the technical point 
of view, both cases can be dealt with in a similar fashion. To study ( p.l| ), we use the 



asymptotic method developed in [g, |], f2(J for the analysis of slow perturbations of one- 
dimensional periodic equations. 



1. The results 

We first recall some of elements of the spectral theory of one-dimensional periodic 
Schrodinger operator Hq; then, we present our results. 



1.1. The periodic operator. For more details and proofs we refer to section |2.1| and 
to [|, 0. 

1.1.1. The spectrum of Hq. The spectrum of the operator Ho defined in ( |0.2|) is a union 
of countably many intervals of the real axis, say [E2 n +i, E 2n+2 ] for n G N , such that 

Ei < E 2 < E 3 < E 4 . . . E 2n < E 2n+ i < E 2n+ 2 < ■ ■ ■ , 
E n — > +00, n — > +00. 

This spectrum is purely absolutely continuous. The points (Ej)j £ ^ are the eigenvalues of 
the self-adjoint operator obtained by considering the differential polynomial ( p.2|) acting 
in L 2 ([0,2]) with periodic boundary conditions (see The intervals [E 2n +i, E 2n+2 ], 

n G N, are the spectral bands, and the intervals (E 2n , E 2n+ i), n G N*, the spectral gaps. 
When E 2n < E 2n+ i, one says that the n-th gap is open; when [E 2n _i, E 2n ] is separated 
from the rest of the spectrum by open gaps, the n-th band is said to be isolated. 

From now on, to simplify the exposition, we suppose that 
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(O): all the gaps of the spectrum of H are open. 



1.1.2. The Bloch quasi-momentum. Let x i— > ip(x,E) be a non trivial solution to the 
periodic Schrodinger equation H ip = Eip such that, for some /x G C and all x G R, 
i/j (x + 1, E) — ^x/; (x, E). This solution is called a Bloch solution to the equation, and // 
is the Floquet multiplier associated to ip. One may write \x = exp(z/c); then, fc is the Bloch 
quasi-momentum of the Bloch solution if). 

The mapping E h-+ fc(.E') is an analytic multi- valued function; its branch points are the 
points Ei, E 2 , E 3 , . . . , E n , .... They are all of "square root" type. 

The dispersion relation k i— > E(fc) is the inverse of the Bloch quasi-momentum. We refer 
to section I2TTI for more details on k. 



1.2. The assumptions on and the analytic continuation of the resolvent. 

We now make assumption (H4) more precise. We introduce the following notation: for 
C > 0, define Cc to be the cone 

C Co = {z G C; C |Imz| < 1 + \Rez\}. 

We assume 

(H4a): there exists Co > such that W : Cc — > C is real analytic, non constant; 
(H4b): there exist (W + , W_) G R 2 and s > 1 such that, for any Ci > C 

(1.1) sup [|z| a |W(s) - W+|] + sup [|z| fl |W(;z) - W-\] < +oo 

Rcz>0 Rez<0 

So IV is short range in cones neighboring +oo and — oo. It is well known that, when Hq 
is the free Laplace operator, this assumption guarantees that the resolvent of H + W 
can be analytically continued from the upper half-plane through the spectrum of Hq as a 



mapping from L 2 Q (R) to L 2 , (M) (for some r > 0) (see [13, and references therein). 
Here, for r G R, we have defined 

l? r = {uG LL(M); e'H«(.) G L 2 (R)}. 

For Imi? > and (el, the resolvent of H^ e at energy E, that is i— > R(E,(,e) : = 
(-ff^e — i?) -1 is a function valued in the bounded operators on L 2 (R); moreover, under 
assumption (H4), a simple resolvent expansion shows that, for any (£"0,(0) G {ImE > 
0} x R, it is analytic in a neighborhood of (£"o, Co)- We prove 

Theorem 1.1. Assume (HI) - (H4) are satisfied. Pick Eq G R sitc/i #ia£ either E — W + G 
o-(#o) \ <M#o), or E — G o-(# ) \ ^(#0), « r ^ 

Then, there exist po > 0, r > 0, > and a complex valued function A : (E,£,e) 1— > 
A(i£, C,e) defined on D(E ,r ) x {(; |/m(| < r } x (0,e ) s,uc ^ /ore G (0,e ), 

• i7ie mapping (E,() \—> A(E,(,e) does not vanish on {E; \E — Eq\ < ro, ImE > 
0} x R; 

• the mapping (E,Q \— > A(E,(,e) is analytic on D(E ,ro) x {(; \Im(\ < r }; 

• the mapping (E,() 1— » A(E,(,e)R(E,(,e) can be continued analytically from 
{E; \E - E \ < r , ImE > 0} x R to {E; \E - E \ < r } x {(; \Im(\ < r } 
as a non-vanishing bounded operator from L 2 o (R) to L 2 po (R); 

• for any E G D(0, ro), the mapping £ 1— > A(i5, £, e) is e-periodic. 
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Following the now classical definition (see [|], [| ^4], pEfl ), we set 

Definition 1.1. Fix ( real and i?o as in Theorem |1.1| . An energy E £ D(Eo,tq) is a 
resonance for if^ e if it is a pole of h-> R(E, (, e) i.e. if it is a zero of E i— >■ A(£", £, e). 

Our goal is to describe the resonances for H^ }E near i?cb an energy satisfying the as- 



sumptions of Theorem Our main interest is in computing the imaginary part of the 
resonances. It is well known that the existence and the positions of resonances will depend 
crucially on the relative position of the spectral window J~{E) := (E — W)(M.) with respect 
to <t(H ) (see @, §). Therefore, one introduces the set W(E) := (E - W)~\(t{H q )). It 
is a subset of the domain Cc where W is analytic. 

1.3. Assumptions on the energy E . Pick E £ R such that either E — W + £ cr(H ) \ 
da(H ), or E — W_ £ cr(H ) \ da(H ), or both hold. A simple Weyl sequence argument 
then shows that E £ a(H^ £ ) (for any ( £ R and e > 0). 

We consider two types of situations and therefore introduce two different assumptions 
(H5): the set W(E ) n R can be decomposed into 

W(E ) ni = U4E ) u ^+(^0) 

where U-(Eq) and U+(Eq) are two-by-two disjoint and they satisfy: 

• U-(Eq) and U + (Eq) are either empty or semi-infinite intervals, respectively 
neighborhoods of — oo and +oo; 

• the finite edges of U-(E ) and U + (E ) (when they exist) are not critical points 
of W. 

and 

(H6): the set W(i?o) H R can be decomposed into 

W(Eq) n r = tZ-^o) u u(E ) u ^+(^0) 

where U-(Eq), U(E ) and U + (E ) are two-by-two disjoint and they satisfy: 

• U{Eq) is a compact interval not reduced to a single point; 

• U^(E ) and U + (E ) are as in (H5). 

• the finite edges of U-(E ), U + (E ) and U(E ) (when they exist) are not 
critical points of W. 

In Fig. H, [l] and ||, we provide some examples of potential profiles; we graphed ( 1— > 
— and, on the vertical axis represented the relevant spectral intervals for Hq. 

Let us now shortly discuss these assumptions. One easily convinces oneself that, if (H5) 
or (H6) holds for some energy Eq, it will hold for all real energies in a neighborhood of 
Eq. Moreover, the fact that the sets U_{E) or U + (E) in the decomposition of W(E) PI R 
are empty will not depend on E (in this neighborhood). Fig. ^provides examples where 
assumption (H5) holds. Under assumption (H6), Fig. ^ provides some examples where 
both U-(E) and U + (E) are non empty, and Fig. ^, examples where either U_(E) or U + (E) 
is empty. The two assumptions (H5) and (H6) correspond to the simplest possible cases: in 
the general case, i.e. when the sets E—W(M.) and ct{Hq) are in a general position, there can 
be more than one compact connected component to the set £ R; E — W(Q £ cr(Ho)}. 



We now state our results on the resonances of H^ £ . 
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(a) (b) 
Figure 1: Examples of W for which assumption (H6) holds and both U± are not empty. 




n+l 




n+1 



(a) 



(b) 



Figure 2: Examples of W for which assumption (H5) holds. 



1.4. Resonance free regions. We begin our results with the energies a neighborhood 
of which do not carry resonances, namely the energies satisfying (H5). Indeed, we prove 



Theorem 1.2 ( fl6fl ). Fix E satisfying (H5). There exist e > 0, 5 > 0, V C C, 
a neighborhood of Eq such that, for e G (0,Eq), for all real (, the set Vq contains no 
resonance of H^^. 

This is quite analogous to what is found in the case when H is the free Laplace operator 
(see e.g. 0,113). 



1.5. Computing the resonances. To study the real energies close to which one finds 
resonances, we introduce a few auxiliary functions necessary to describe the resonances. 
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(a) 



(b) 



Figure 3: Examples for which assumption (H6) holds and either C/_ or U + is empty 



1.5.1. The complex momentum and its branch points. The complex momentum ( i— > k(Q 
is defined by 

(1.2) K(() = K((,E) = k(E-W(0). 

As k, k is analytic and multi-valued. As the branch points of k are the points (Ei) ie ^, 
the branch points of n satisfy 

(1.3) E-W(£) = E j , j GN*. 

As E is real, the set of these points is symmetric with respect to the real axis, to the 



imaginary axis and it is 27r-periodic in (. More details are given in section 2.2.1 



Fix Eq satisfying assumption (H6). By ( |1.3| ), for E real close to Eq, the ends of U(E) 
and the finite ends of U_(E) and U + (E) (when they are not empty) are branch points of 
k. To fix ideas, let us define 

U(E) = {C (E),ti(E)l U.(E) = (-oo, C(E)] and U+(E) = [C + (E], +oo) 

When U + (E) (resp. U-(E)) is empty, we set C + (-^) = +°° (resp. (~(E) = —oo). 
One shows that there exists a determination of the complex momentum, say k,q, and a 
real neighborhood of E , say V^, such that, for E e V R , « (C) ^ [0, n] for ( e U{E). One 
defines 

rti(E) 1 
(1.4) * (E)= ^(C,E)d( and 6k = -[k((+(E), E) - k (Q(E), E)}. 

One shows 



Lemma 1.1 ([IB|). The constant 5k belongs to { — 1,1,0} and is independent of E in 
V R . The function $o : Vq R — > ^ + can be extended to a function analytic in a complex 
neighborhood of Eq. Moreover, there exists C > such that, in this neighborhood, one has 



;i.5) 



& (E)\>- and & (E)-6k>0. 



Assume U-(E) or U + (E) or both are not empty. One then proves that the imaginary 
part of the determination k keeps a fixed sign in the interval between U_(E) and U(E), 
and between U(E) and U+(E), i.e. on the intervals ((-{E), (o(E)) and ((£(E),(+{ E ))- 
We define 

r(o( E ) ( rC+W > 

(1.6) S-(E) = ± Im/<CR resp. S+(E) = ± Imre(C)dC 

where ± is chosen so that S+(E) and S-(E) are positive. One proves 

Lemma 1.2 ( ||16|| ). Pick v G {+,—}■ The function S v : Iq R — > K + can 6e extended to a 
function analytic in a complex neighborhood of E . 

When U + (E ) (resp. U-(E )) is empty, it stays empty for E close to E , and one sets 
S + (E) = +00 (resp. S-(E) = +00). Define the tunneling coefficients 

(1.7) t±{E) = e - s±{E)/£ and t(E) = t + (E) + t_(E). 
1.5.2. Resonances. We prove 

Theorem 1.3 (|T§). Fix E satisfying (H6). There exist e > 0, 5 > 0, V C C, 
a neighborhood of Eq, and a real analytic function E > 3>(E,e) defined on Vq x [0, £0] 
satisfying the uniform asymptotics 

(1.8) ${E,e) = $(E) + 0(e) when e -> 0, 

sitc/i t/iai, i/ one defines the finite sequence of points in Vq H R, sa?/ (-E'Oi := (E l ((,e))i, 
by 

(1.9) e) = -tt5kC+ - +tt/, ZeZ, 

i/ien, /or e G (0,£ ), /or a// real (, the resonances of Hq £ in Vq are contained in the union 
of the disks 

(1.10) D l := {E G V Q ] \E - E l \ < e~ 5o/e } . 

More precisely, each of these disks contains exactly one simple resonance, say E l ((,E), 
that satisfies 

(1.11) Im(E l ((, E)) = e ■ Cq ■ t(E l )(l + o(l)) 
where o(l) — > when s — * and Cq is a constant. 



Let us discuss the resonances obtained in Theorem |1.3| , in particular, their behavior 
as functions of £. By Theorem |1 . 1| , the function ( > A(E,(,z) is £-periodic and by 
Theorem |1.3| , the resonances are simple, the resonances being the zeroes of E 1— > A(E, (,e), 
they are £-periodic functions of (. 

Theorem |1.3| says that their imaginary part does not oscillate to first order. It also 
gives some information about the oscillations of the real part. Therefore one has to 
distinguish the two cases 5n 7^ and 5k = 0. 

When 5k 7^ 0, the quantization condition ( |1.9| ) and equations (|1.5|) and ( |1.9|) show 
that the real part of the resonances is monotonous in Q. But, as they come in a sequence, 
renumbering the resonances, they can also be seen as oscillating with amplitude roughly 
Ce (for C > 0). Of course, to see this oscillation phenomenon, one needs to renumber the 
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resonances. The situation is similar to that encountered when studying Stark- Wannier 
resonances (see in this case also, the resonances still exhibit oscillations of frequency 
27r/e but the amplitude is e (the Stark- Wannier ladders). 

When 5k, ^ 0, Theorem |1.3| shows that the oscillations of the real part are at most 
exponentially small. Under more precise assumption than those made in the present talk, 
one can compute the amplitude of the oscillations of the resonances ([0). An analogous 
behavior was found for the eigenvalues of Hq £ outside the essential spectrum in |2t|; they 
oscillate with a frequency 2ir/s and an amplitude that is exponentially small in e; the 
amplitude was given by a complex tunneling coefficient. 

Moving the energy E (for a fixed W), one sees that, in some cases, one can pass 
continuously from assumption (H5) to (H6) and vice-versa. It would be interesting to 



study what happens to the resonances uncovered in Theorem [L3] when one crosses this 
transition. Such a study has been done in O, [13] in the case when V = 0. 



1.6. The heuristics. We now discuss the heuristics explaining these results. In the 
figures below, in the part indexed a), we represented the local configurations of the po- 
tential to show the assumption made on the energy E in the different cases. In the part 
indexed b), we represented the phase space picture of the iso-energy curve i.e. the sets 
{(k, C) G K 2 ; k = k(E — W(())}. It is 27r-periodic in the ^-direction (i.e. the vertical 
direction) so we depicted only two periods. On the same picture, we also show some loops 










(a) 



(b) 



Figure 4: The phase space picture 



that are dashed. These are special loops in {(k, () G C 2 ; k = k(E — W(())} that join the 
connected components of {(k, Q G M 2 ; k = k(E — W(Q)}. 

Assumption (H5) just means that the real iso-energy curve is empty, in which case there 
are no resonances as we saw. We now assume (H6) holds. 

In Fig. ||, the compact connected components of the iso-energy curve consist of a 
single torus per periodicity cell. The torus gives rise to the quantization condition (|1.8|) : 
the corresponding phase is obtained (to first order) by integrating the canonical one- 
form nd( along this torus. The fact that this corresponds to a resonance rather than 
to a eigenvalue is due to the fact that the iso-energy curve has components leading to 
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infinity. The dashed curve represents the instanton linking the two. To compute the 
action determining the lifetime of the resonance, one integrates the canonical one-form 
ndC, along this curve. 




(a) (b) 



Figure 5: The phase space picture 



In Fig. the compact connected components of the iso-energy curve are absent. 
Nevertheless, one can compute a phase (to first order) by integrating the canonical one- 
form ndC, along a period of the central curve i.e. the one that is extended in the K-direction. 
When the connected component extended in the ^-directions are absent, one can show 
that such a phase does not give rise to eigenvalues. The presence of these components 
leading to infinity in the ^-direction gives rise to resonances. Again, to compute the action 
determining the lifetime of the resonance, one integrates the canonical one-form ndC, along 
the dashed curve. 



2. The adiabatic complex WKB method and the proof of Theorem 1.1 



In and ||, we have developed a new asymptotic method to study solutions to 
an adiabatically perturbed periodic Schrodinger equation i.e., to study solutions of the 
equation 

(2.1) C) + + Wi - £X + C)) ^ (x ' C) = E ^ C) 

in the limit e — ► 0. The function ( i— > W(() is an analytic function in a neighborhood of 
the real axis. The main idea of the method is to get the information on the behavior of the 
solutions in x from the study of their behavior on the complex plane of (. The natural 
condition allowing to relate the behavior in ( to the behavior in x is the consistency 
condition 

(2.2) t/;(x + l,(,E)=iP(x,( + e,E). 

One can construct solutions to both ( |2.1| ) and ( |2.2| ) that have a simple asymptotic behavior 
on certain domains of the complex plane of £. 

The ideas underlying this result are simple. For solutions of ( |2.1| ) satisfying ( |2.2| ), it is 
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equivalent to know their behavior for large x or for large (. Hence, the idea is to study 
their behavior for large (. Keeping x fixed to some compact interval, say [—1, 1], to first 
order in e, a solution to (12. II) is a solution to 



d 2 

(2.3) -—ip(x) + V(x)ip(x)=£ip(x), x6l, 

for £ = E — W((). Hence, the Bloch solutions are special solutions to this equation. From 
these solution, multiplying them by a suitably chosen function of £ (i.e. independent of 
x), one can construct solution to both ( ^.1| ) and Q2.2| ) at the same time. The standard 
form of the solutions is given below in ( [2.11| ). 



We first describe the standard asymptotic behavior of consistent solutions; in order to 
do this, we recall some fact on periodic Schrodinger operators on R. Then, we review the 
theory developed in || |J on the existence of consistent solutions with these asymptotics. 
We next describe some new results used to control the behavior of the solutions at infinity. 



We conclude this section with the proof of Theorem LI . 



2.1. Periodic Schrodinger operators. In this section, we discuss the periodic Schrodin- 
ger operator ( p.2| ) where V is a 1-periodic, real valued, L^ oc -function. First, we collect 
well known results needed in the present paper (see O, |^, [Tj| |2T], [23|]). In the second part 



of the section, we introduce a meromorphic differential defined on the Riemann surface 
associated to the periodic operator. This object plays an important role for the adiabatic 
constructions (see 131). 



2.1.1. Bloch solutions. Let ip be a nontrivial solution of the equation ( |2.3| ) satisfying the 
relation ip (x + 1) = X (x) for all x G R with A G C independent of x. Such a solution is 
called a Bloch solution, and the number A is called the Floquet multiplier. We now discuss 
properties of Bloch solutions (see |p"Tf ). 



As in section we denote the spectral bands of the periodic Schrodinger equation by 
[Ei, E 2 ], [E 3 , E4], . . . , [E 2n +i, -EW1+2], .... Consider S±, two copies of the complex plane 
£ G C cut along the spectral bands. Paste them together to get a Riemann surface 
with square root branch points. We denote this Riemann surface by S. In the sequel, 
7r c : S 1— > C is the canonical projection. 

One can construct a Bloch solution ip(x, £) of equation ( |2.3j ) meromorphic on S. For any 
£, we normalize it by the condition ifj(l,£) = 1. Then, the poles of £ 1— > ip(x,£) are 
projected by ir c either in the open spectral gaps or at their ends. More precisely, there 
is exactly one simple pole per open gap. The position of the pole is independent of x 
(see 0). 

Let : : S 1— ► S be the canonical transposition mapping: for any point £ G S, the point 
£ is the unique solution to the equation ir c {£) = E different from £ outside the branch 
points. 

The function x t— > tp(x,£) is one more Bloch solution of ( ^.3| ). Except at the edges of 
the spectrum (i.e. the branch points of S), the functions ip(-,£) and i/j(-,£) are linearly 
independent solutions of (|2.3|). In the spectral gaps, they are real valued functions of x, 



and, on the spectral bands, they differ only by the complex conjugation (see [IT 
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2.1.2. The Block quasi-momentum. Consider the Bloch solution ip(x, £). The correspond- 
ing Floquet multiplier A (£) is analytic on S. Represent it in the form X(£) = exp(ik(£)) . 
The function k(£) is the Bloch quasi-momentum. 

The Bloch quasi-momentum is an analytic multi- valued function of £. It has the same 
branch points as ijj{x,£) (see |TT]j ). 



Let D G C be a simply connected domain containing no branch point of the Bloch quasi- 
momentum k. On D, fix ko, a continuous (hence, analytic) branch of k. All other branches 
of k that are continuous on D are then given by the formula 

k ±t i(£ ) = ±k (£) + 2ttI, I G Z. 

All the branch points of the Bloch quasi-momentum are of square root type: let E\ be 
a branch point, then, in a sufficiently small neighborhood of Ei, the quasi-momentum is 
analytic as a function of the variable \f£ — E~f, for any analytic branch of k, one has 



k(£) = k l + c l y/E-E l + 0(£-Ei), ci^O, 
with constants fc; and q depending on the branch. 

Let C + be the upper complex half-plane. There exists k p , an analytic branch of k that 
conformally maps C + onto the quadrant {k G C; Imk > 0, Rek > 0} cut along compact 
vertical intervals, say irl + il\ where I G N* and J; C R, (see ||11|| ). The branch k p is 
continuous up to the real line. It is real and increasing along the spectrum of Hq; it maps 
the spectral band [E 2 n-i, E 2n \ on the interval [n{n — l),7m]. On the open gaps, Rek p 
is constant, and lmk p is positive and has exactly one maximum; this maximum is non 
degenerate. 

We call k p the main branch of the Bloch quasi-momentum. 

Finally, we note that the main branch can be analytically continued on the complex plane 
cut only along the spectral gaps of the periodic operator. 

2.1.3. Meromorphic differential Q. On the Riemann surface S, consider the function 

JqiI>(x,£) (ip(x,£) — ik(£)xvjj(x,£)) dx 



(2.4) u{£) - 

J ijj(x,£)i/j(x 1 £)dx 

where k is the Bloch quasi-momentum of ip, and the dot denotes the partial derivative 
with respect to £. This function was introduced in |10j (the definition given in that paper 



is equivalent to (|2.4j )). In |fL0|j , we have proved that u has the following properties: 

(1) the differential Q = ud£ is meromorphic on S; its poles are the points of P U Q, 
where P is the set of poles of ip(x, £), and Q is the set of points where k'{£) = 0; 

(2) all the poles of Q are simple; 

(3) Vp G P\Q, res p Q = 1; G Q \ P, res q Q = -1/2; WePnQ, res r Q = 1/2. 

(4) if ix c (£) belongs to a gap, then ui{£) G 1R; 

(5) if n c (£) belongs to a band, then uj{£) = uo{£). 

2.2. Standard behavior of consistent solutions. We now discuss in more detail two 
analytic objects central to the complex WKB method, the complex momentum defined 
in ( |1.2[ ) and the canonical Bloch solutions defined below. For ( G T>(W), the domain of 
analyticity of the function W, we define 

(2.5) £(Q=E-W(C) 

li 



The complex momentum and the canonical Bloch solutions are the Bloch quasi-momentum 
and particular Bloch solutions of the equation 

(2.6) _A_^ + v^ = £( C ty. 
considered as functions of C- 

2.2.1. The complex momentum. For C G T>(W), the domain of analyticity of the function 
W, the complex momentum is given by the formula k(Q = k(S(Q) where k is the Bloch 
quasi-momentum of ( P-2| ). Clearly, k is a multi- valued analytic function; a point ( such 
that W'(C) 7^ is a branch point of k if and only if it satisfies ( |1.3| ). All the branch points 
of the complex momentum are of square root type. 

A simply connected set D C T>(W) containing no branch points of k is called regular. Let 
k p be a branch of the complex momentum analytic in a regular domain D. All the other 
branches analytic in D are described by 

(2.7) = ±k p + 2irm where m G Z. 

2.2.2. Canonical Bloch solutions. To describe the asymptotic formulae of the complex 
WKB method, one needs Bloch solutions of equation (|2.6|) analytic in C on a given regular 



domain. We build them using the 1-form Q = to d£ introduced in section 2.1.3 



Pick Co, a regular point i.e. a point that is not a branch point of k. Let Sq = S(Q 



Assume that £ G" P U Q (the sets P and Q are defined in section |2.1.3| ). In Uq, a 
sufficiently small neighborhood of Sq, we fix k, a branch of the Bloch quasi-momentum, 
and ip±(x, S), two branches of the Bloch solution ip(x, S) such that k is the Bloch quasi- 
momentum of Also, in Uq, consider Q±, the two corresponding branches of Q, and 
fix a branch of the function £ t— > q{£) = ^k'(S). Assume that Vq is a neighborhood of Co 
such that £{Vq) C Uq. For ( G Vo, we let 

(2.8) *±(x, C) = q(S) e J "o n± <0±(x, £), where £ = 8(C). 

The functions \I/± are the canonical Bloch solutions normalized at the point Co- The 
quasi-momentum associated to these solutions is = k(E — W(()). 

The properties of the differential Q imply that the solutions ty± can be analytically con- 
tinued from Vq to any regular domain D containing Vq. 
One has (see 0) 

(2.9) w(^ + (; C), *_(-, 0) = Co), *-(•, Co)) = fc^M^+O, £>), ^-(-, £o)) 
As £ ^ <3 U {-E;, Z > 1}, the Wronskian w($+(-, (), #_(•, C)) does not vanish. 

2.2.3. Solutions having standard asymptotic behavior. Fix E = E . Let D be a regular 
domain. Fix Co G -D so that £(Co) ^ PUQ. Let k be a branch of the complex momentum 
continuous in D, and let \&± be the canonical Bloch solutions defined on D, normalized 
at Co and indexed so that k be the quasi-momentum for 
We recall the following basic definition from || 

Definition 2.1. Fix r\ G {+, — }, £q E D, X > 1, Vq, a complex neighborhood of some 
energy E and D a domain in the C-plane. We say that /, a solution of ( |2.1|) , has standard 

asymptotics (or standard behavior) / ~ exp(r)^ ndu) ■ ^ v in (— X, X) x D x Vq if 
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• / is denned and satisfies fl2TT| ) and (|2.2|) for any (x, C, -E) € (—X, X) x D x V ; 

• / is analytic in £ G -D and in i£ G Vo; 

• for any K, compact subset of D, there exists V C Vo, a neighborhood of £b, such 
that, for (x,(,E) G (— X, X) x K x V, f has the uniform asymptotics 

/(*, C, E, e)=e V " & E) du (%(x, £ E) + g v (x, £ £, e)) 

where 

lim sup |^(a;, C, -E 1 , e)| = 0. 
xe(-x,x) 

• this asymptotic can be differentiated once in x. 

Let (/_)_, /_) be two solutions of ( F2.1|) having standard behavior f± ~ e £ J * in 
(— X, X) x D x Vq. One computes 

w(f + , /_)=«;(*+,*_) + o(l). 

By ( p.9| ), for ^ in any fixed compact subset of .D and £ sufficiently small, the solutions 
(/+,/_) are linearly independent. 

2.3. The main theorem of the adiabatic complex WKB method. A basic and 
important example of a domain where one can construct a solution with standard asymp- 
totic behavior is a canonical domain. Let us define canonical domains and formulate the 
basic result of the adiabatic complex WKB method. 

2.3.1. Canonical lines. We recall that a curve 7 is vertical if it intersects the lines {Im ( = 
Const} at non-zero angles 6, < 9 < ir. Vertical lines are naturally parametrized by 
ImC- 

\ set D in the ^-plane is said to be bounded regular at energy Eq if D, the closure of D, 
is compact and does not contain any branch point of C 1 — >■ k(£, Eq). 

Let 7 be a C l bounded regular vertical curve at energy Eq. On 7, fix a point Co an d k, a 
continuous branch of the complex momentum. 

Definition 2.2. The curve 7 is canonical with respect to k at energy E if, along 7, for 
y — Im£, one has 

— lm([ k(u, E )du) > and — (im / (k(u, Eq) — 7r)du ) < 
dy \J( J dy V ico / 

One easily checks that, 

• if a vertical curve is bounded regular at energy Eq, it is bounded regular at any 
energy E in a neighborhood of Eq; 

• if a vertical curve is bounded regular and canonical with respect to k at energy 
Eq, it is also canonical with respect to k at any energy E in a neighborhood of Eq. 

13 



2.3.2. Canonical domains. Let K be a bounded regular domain at energy Eq i.e. K is 
compact and does not contain any branch point of C i— > So). On K, fix a continuous 
branch of the complex momentum, say k. 

Definition 2.3. The domain K is called canonical (with respect to k at energy Eq) if 
there exists two points £1 and C2 located on the boundary of K, such that .fT is the 
union of curves connecting Ci and C2 that are canonical (with respect to k at energy Eq). 

2.3.3. T/ie main theorem of the adiabatic complex WKB theory. One has 

Theorem 2.1 (|6|, |8|]). Lei K be a bounded domain canonical with respect to n at energy 
Eq. Fix X > 1 and Co £ Then, there exists Vq, a neighborhood of Eq such that, for 
sufficiently small positive e, there exist (f±), two solutions of ( p.l|) , having the standard 
behavior in (—X,X) x K x Vq that is 



C 
Co 



For any fixed £/ie functions (C, E) 1— > /±(x, £, i?) are analytic in S(K) x Vo where 

S(K) is the smallest horizontal strip containing K. 

In addition to this result, we need analogues for domains that are neighborhoods of 



infinity. Such results were developed in |19], |20| ; they are not sufficient for the purpose of 
the present paper. Hence, we present two new results in the next section. 

2.4. Consistent solutions at infinity. The standard behavior at infinity is given by 



Theorem 2.2 (|16[). Assume (H1)-(H4) are satisfied. Fix rj £ {+, — } and a £ (0,1). 
Pick Eq £ E such that either Eq — W^E (t(Hq) \ da(H ). Fix X > 1 and C\ > C (where 
Cq is defined in assumption (H4), section Then, there exist C > and 5 > such 

that, for any E £ D(E ,5), the cone C T] = {( £ C; r]Re( > C, Ci\Im(\ < r]Re(} is 
a canonical domain. More precisely, fix ( £ C T] fl K. Then, there exist a branch of the 
quasi-momentum n(Q = E) that satisfies 

(2.10) V(C, E) £ C r] x D(E , S), < Re E) < tt 

and a function (x, (, E) 1— > f(x, (, E) defined on (—X, X) x C v x D(E , 8) that satisfies: 

• it is a consistent solution to equation ( |2.1| ) 

• for x £ (— X, X), ((, E) 1 — > f(x, C, E) is analytic on C v x D(E , 5); 

• the function f admits the asymptotics 

(2.11) /(*, C, E) = e^ J co < u ^ du (Vfa E) + g(x, C, E, e)) 

where ^ is the canonical Bloch solutions associated to k (see ( |2.8|) ) and 

(2.12) sup sup I^CM^C^e)! < +00. 

ee(0,£ ) CGC+ 

E&D(E ,5) 

• the asymptotics can be differentiated once in x; 



• define the function (x, £, E) 1— > f*(x, (, E) = f(x, (, E) ; then, (/, /*) form a basis 
of consistent solutions and satisfy 

(2.13) W(f(; C, E), f*(; C, E)) = w(V + (; Co), Co))- 

14 



We will use a second result concerning consistent solutions near oo. We won't need a 
basis in this case: a single solution (actually the Jost solution) will be sufficient. It is 
given by 

Theorem 2.3 (fL6|). Assume (H1)-(H4) are satisfied. Fix 77 G {+, — } and a G (0,1). 
Pick E G M such that E — W v G" a(H ). Fix X > 1 and C\ > Cq (where Cq is defined 
in assumption (H4), section \TJ\ ). Then, there exist C > and 5 > such that, for 
any E G D(E ,5), the cone C v = {( G C; r]Re( > C, Ci\Im(\ < r]Re(} is a canonical 
domain. More precisely, fix(o G C^flM. Then, there exist a branch of the quasi-momentum 
K (C) — K ((, E) that satisfies, for some a G {0, 1} ; 

(2.14) V(C-E) G C v x D(E Q ,S), 0<ImK((,E) and k*((, E) = 2ixa — k((, E). 

and a function (x, (, E) 1— > f(x, (, E) defined on (—X, X) x C v x D(E , 5) that satisfies: 

• / is a consistent solution to equation (|2.1| ) 

• for x G (— X, X), {C,,E) 1 — > f(x,(,E) is analytic on C v x D(E ,5); 

• the function f admits the asymptotics 

(2.15) f(x, C, E) = e"^c C «(".*>*• C, E) + ^(x, C, £?, e)) 

where (\l/ + ,\I / _) are i/ie canonical Bloch solutions associated to k (see (|2.8| ) j and 



• ^/ie asymptotics can be differentiated once in x; 

• the functions satisfy 

(2.16) V(x,C,S) G (-X,X) x x D(E ,5), f\x,t,E) = e~ 2 ^-^f(x, (, E). 



In PLQfl ■> a result analogous to Theorem is proved, the main difference being the 
region in ( where the statements hold. 

2.5. The continuation diagrams. Now to achieve the goal described in section 0, 
namely, to construct consistent solutions to ( |2.1| ) with known asymptotics in large com- 
plex domains of £, as we have the form of standard asymptotics locally both at a finite 
point in ( and near infinity, we only need patch together the various domain on which we 
found consistent basis with standard asymptotics. An alternative way to obtain global 
asymptotics was developed in || |TI], |7|, [J. It consists in the proof of quite general contin- 
uation lemmas that describe geometric situations in the complex plane of £ under which 
one can prove that, starting from a canonical domain and solutions with standard asymp- 
totics in this domain, one can continue them. The main obstacles to continuation (i.e. to 
the validity of standard asymptotics) are the nodal lines of Im k and the branching points 
of k. We will not discuss this further and refer to [16| for the details. 
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